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Abstract 

^ , It is shown that the interaction of hehcity-1 waves of gravity and matter in 

^ \ a thin slab configuration produces new types of instabihties. Indeed, a spin- 

2 hehcity-l mode interacts strongly with the shear motion of matter. This 
mode is unstable above a critical wavelength, Ac = -sJ-kc^ /2Gp. This should 
be compared with Jeans wavelength, Aj = y^vrcf/Gp, where Cg is the sound 
speed. The two instabilities are of course different. For the case analyzed, a 
plane parallel configuration, Jeans instability appears through a density wave 
perturbation, the material collapsing into a set of plane-parallel slabs. On the 
other hand, the helicity-1 wave instability induces a transverse motion in the 
fluid that tends to shear in the material along the node of the perturbation. 
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Acoustic waves are affected by gravitational forces. One of the consequences is 
Jeans instability [Q (see also, e.g., [0, produced above a certain critical wave- 
length. There is an electromagnetic counterpart in plasma density waves. In the 
hmit of large wavelengths their frequency square tends to a finite positive value, and 
the plasma acquires the plasma frequency, whereas in the Jeans case the frequency 
square turns negative. This is noted by Weinberg |Q but he does not go further. 
Plasma density waves can be seen as longitudinal electromagnetic waves which are 
not possible in vacuum but are present in media with the plasma frequency hav- 
ing the role of a photon mass. By the same token acoustic waves can be identified 
with longitudinal waves of the gravity field but now, instead, they possess an imag- 
inary mass. Transverse electromagnetic waves, in addition to longitudinal waves, 
are also modified in a plasma and acquire a mass equal to the plasma frequency. 
The questions posed in this paper concern waves of the gravity field, or wave metric 
perturbations, propagating in matter. More specifically, we are interested in the 
following question: are there waves which get a finite real mass or do they get an 
imaginary mass leading to instabilities other than Jeans? 

Waves of the gravity field can be classified according to their spin and helicity. 
Here we are interested in spin 2 waves which can have helicities 2, 1 and 0. Helicity- 
2 waves represent radiative gravitational waves and are simply called gravitational 
waves. Helicity-1 waves, as well as helicity-0, do not exist in vacuum but can prop- 
agate in matter. The propagation of waves of the gravity field through matter has 
been studied in two main contexts, namely, cosmological and astrophysical. Within 
the cosmological context these waves appear as perturbations of the background 
universe (see also 0) and there is now a growing interest in the primordial grav- 
itational wave background since it seems possible to detect it in the future ||^, |[. 
Within the astrophysical context they have been studied mainly in relation to stars 
(see e.g., [^). Usually one considers these stars as either spherically or axisymmetric. 
The waves in the star, i.e., perturbed quadrupolar or higher polar modes, appear as 
convective velocity perturbations and density perturbations which through the stel- 



lar boundary match with gravitational waves propagating in vacuum |T^, Still 
within the astrophysical context, there has been another line of research, which is 
the one that most concerns us here, analyzing non-dissipative |P, 0, |l^ and dissipa- 
tive propagation of waves through matter |]14|, |15|, [l^, H . The results obtained so 
far confine themselves either to the short wavelength limit or to transverse helicity-2 



gravitational wave modes (e.g., |T^, |T6|)- In either case the interaction of wave and 
matter is minute, and usually the wave propagates undisturbed. 

In this paper we want to go a step further, and reach the level where the wave- 
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length, A, is of the order of or larger than the background curvature, ~ c/ \fGp. This 
is not possible to achieve for spherical or quasi-spherical geometries, since it implies 
systems with a size of the order of their Schwarzschild radius. However, extended 
finite thin plates either disk-like or parallelepiped-like are amenable to the large A 
regime. In order to progress we will study a specific case, namely a wave propa- 
gating through a plane parallel slab. We are then able to find some modes which 
propagate and truly interact with matter. The transverse helicity-2 gravitational 
wave mode propagates freely without interaction as has been observed in several 
papers (see e.g. |1^). On the other hand, a transverse helicity-1 shear wave 
mode interacts strongly with the shear motion of matter. Indeed, above a certain 
critical wavelength this mode is unstable, as we will now show. 

Consider a fluid with proper energy density p, pressure p and with a given equa- 
tion of state Consider then it has suffered a small departure from a state where 
it is at rest in a Minkowski space with homogeneous density, po? and pressure, po- 
The fluid energy-momentum tensor is (c = 1) 

T^,u = ip + p)u^u^-pg^u , (1) 

where denotes the 4- velocity vector. Greek indices /x, u, etc., take values (or 
t), and i, with i = 1,2,3 (or x, y, z, respectively). The metric g^^i, is denoted in 
the Minkowski limit as r/^i, = diag(l, —1, —1, —1). The fluid and the gravitational 
field form a coupled system. The change of T^^, due to a fluctuation in the metric is 
determined by the conservation laws of energy and momentum, V^T^y = and in 
turn, the gravitational field is determined by the Einstein equations 

R^, = SnGiT^, - \g^,T) . (2) 

One has to distinguish between a static gravitational field, g^^u = Vfiu + ^fj,u, caused 
if not for anything else by the fiuid at rest, from the wave modes /i^i, defined as 

Qfiu = ri^y + h^y + h^y . (3) 

The field h^i, is related through Einstein equations, 

R^4h] = SnGif^, - ^g^^T) , (4) 

to the energy- momentum tensor, T^^u, of the fiuid in static equilibrium under the 
external field hf^u itself. T^^u is characterized by p and p, which in turn deviate 
slightly from po and po, see below. 
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We adopt the de Bonder gauge (see, e.g., |T8|, |T9|, ^ 



5^ V - = , (5) 



for both h^jj and h^y. Here and in the following the Minkowski metric rj^^ is used 
to raise and lower indices of hfj,u and d^. In this gauge the equations of energy- 
momentum conservation are in leading order in the gravitational field 

The limit h^^^ = gives the conservation equations of T^j^y. After subtracting that 
static contribution one obtains the equations governing the variations t^^ = 5T^y = 
T^y — T^y induced by the fluctuations h^y. In first order 

d'r,y - dyh^'^ = , (7) 

where Tj^J = diag{po, Po, Po, Po) replaces T^j,, an approximation that will be made 
whenever T^^y is multiplied by the field hajj because T^y — Tj^J is already suppressed 
by hap ■ Equations (0) give the linearized equations (where, '= d/dt), 

6p - (po + Po) Ui^i + ^(po + Po) hoo = , (8) 

(Po + Po) Ui - 5p^i - ^(po + Po) hoo,i = , (9) 

which yield r^j, as linear functions of ha/B when considering periodic waves. This 
system is completed by Einstein equations (|]) and (^ resulting in the approximate 
equality, 

R,y[h + h]- R,y[h] = 8nG (r^y - \r^^y{r - T^^^ h^") - '^T^'^h^y) , (10) 

where the traces are defined as r = r^yrj^^'^ and T'^'^^ = T^S) r]^^ . In the right-hand 
side we neglect terms that go as, omiting the indices, Thh^ {T—T^^^)h or r/i, because 
they are suppressed by extra factors of h^j_y and henceforth only appear at order. 

In the weak field limit the left-hand side of Eq. (|1^) is approximated by —^Ohfj_y 
whereas in the right hand side the explicit dependence of the matter terms in hfj_y 
is neglected. In the study of the propagation of gravitational waves through matter 
one usually considers the high frequency regime, where the wavenumber k is much 
greater than the background curvature, i.e., k"^ ^ SnGp, enabling one to neglect 
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in the equations the matter density p. However, matter effects become important 
when k"^ approaches the background curvature, i.e., fc^ ~ SirGp. For such small k 
higher order terms of the Ricci tensor, such as A/ioo hf^i, ~ AnGp h^^, are of the same 
order as k^h^y and still linear in h^y. This calls for a higher order expansion of the 
Ricci tensor in equation ( |TOD up to terms linear in both h^^ and background field 
h^y. Such an expansion includes products of h^y derivatives and h^y derivatives and 
products of h^i, times second derivatives of h^i, so that translational invariance is 
lost and monochromatic waves cease to be in general solutions of Einstein equations. 
This may be circumvented in some cases, such as in the solutions presented below. 

Consider a thin wall of fiuid lying over the plane a; = with translational symme- 
try along the y — z plane. If the pressure is everywhere much smaller than the density 
the Einstein equations @} admit in first order a Newtonian solution h^y = 20(x) 6^^, 
with 0" = AirGp. Consider in addition that the fiuid is in a nearly incompressible 
state, i.e., with a practically uniform density p ~ po- A solution with refiection 
symmetry about the plane x = possess an acceleration 0' = AnGpox and pressure 
p = 27iGpl{£'^ /A — x^), p vanishing at the boundary planes x = ±i/2. The density 
variation is negligible across the wall if Ap <^ c^po for a given speed of sound c^. 
That is true if the wall thickness i obeys £ <^ Xj, where Xj = y^iTcl/Gp is the Jeans 
wavelength. Note that the non-relativistic condition p <^ p is a weaker condition 
since it implies £ -C 1/ y/Gp. The magnitude of the potential (j) at the wall depends 
upon the longitudinal dimension of the wall, L, because at distances larger than L 
it decays with the inverse distance to the wall. One obtains |0| ~ GpoiL at the 
wall. Therefore the Newtonian condition yields iL <^ 1/Gp. To summarize, 

the background we are assuming consists of a fiuid with uniform density and small 
pressure associated with a Newtonian gravity potential. Other details, such as the 
precise fall off of the pressure at the rim of the wall, are no needed here. 

To obtain from Eq. (|T0|) an accurate system of linear equations in /i^j^ up to 
order of Gp it is sufficient to have in leading order of Gp in the left-hand 
side of Eqs. ([T0|) . There are terms proportional to 4>'h^y^a and (ph^y^ais- For 

definiteness we consider waves propagating along the z axis, i.e., r^^ and h^y varying 
as exp(— itut + ikz). It is conceivable that for large enough wavelengths the terms 
proportional to (f)' and cj) are much smaller than the terms with 0" = AirGpo so 
that solutions of that type may be found by neglecting the 0' and terms. It is not 
enough however that /c^0 and k(f)' be much smaller than 0". The three kinds of terms 
mix in the Bianchi identities and therefore one should not expect that the Einstein 
equations remain always consistent with each other if some terms are neglected. 
Nevertheless two modes were found that provide consistent solutions of Einstein 
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equations with definite wavelength along the z axis. They are hyi and /102 — ^ 



23- 



After neglecting the 0' and terms (and pressure) Eqs. (10) give for the h\2 
mode 

(^2 _ ^2 _ 2<\>")h^2 = 167rG'(ri2 - ^Pohu) , (11) 

while the fluid equations of motion (§) and (P) yield r^u = 0. Hence, the hi2 mode 
satisfies u;^ = k"^ and propagates at the vacuum speed of light without disturbing 
the fluid. This absence of dispersion is in agreement with various previous results 



regarding the propagation of transverse gravitational waves in media (see, e.g., |T4 



For the /102 and /i23 mode (A;/i23 = — t^/io2 in the gauge (|^)) Eqs. (p!0[) give, after 
neglecting the 0' and terms as before: 
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[UJ 



2 U2 



e)ho2 = -167rG(5T°2 - lpo/io2) , (12) 
P)/i23 = 167rG(5r23 - ipo/i23) , (13) 



where ^T"^ = po^^ is the momentum density along the y— direction and the fluid 
velocity is = — M2 + /io2- The terms proportional to po^o2 and po^23 are due 
to the background curvature Gpo- The fluid equations of motion (||) and yield 
Ui = = 6p = 6p hence, = /io2 and ST"^^ = 0. The resulting dispersion relation is 

= fc2 - SnGpo . (14) 

In contrast to the hu, this /io2 — ^23 mode exhibits an instability at very large 
wavelengths, k'^ < SirGpo, that is induced by matter. Most notably, the polarization 
carries the helicity states ±1 and not ±2, in contrast to the gravitational waves that 
propagate in vacuum. In vacuum, of course, this mode can be gauged away through 
a coordinate transformation but that is no longer true in matter as can be seen by 
the following argument: /i23 can be eliminated with a transformation of the type 
y y + ^{t,z) but /io2 transforms into /iq2 = (1 ~ ^^/^^)^02, which is not zero, 
because the phase velocity uj/k is different from unit in matter. Alternatively, one 
may gauge away u"^ and ho2 but not h23. It is easy to prove that no coordinate 
transformation can bring simultaneously all the components h^j^^, and -u* into zero. 
The motion of the fluid can be characterized in a gauge independent way by the 



covariant derivative of the velocity, u^-i^ (see |£T|). In leading order, the important 
non-zero covariant derivatives are 

U2;3 = M3;2 = "^(^3^02 " <9o/?.23) • (15) 
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It means that this is a shear motion, u^-^y = a^^^ with the square of the shear 
magnitude defined as the gauge invariant quantity cr^ = a^uO^^ l2. In the gauge we 
use, 

(y = U2;3 = p O3ho2 , (16) 

which shows again the crucial role played by matter. These helicity-1 waves can then 
be called shear waves. The components m^, /io2, ^23 correspond to what are called 
vector perturbations in cosmological perturbation theory They are interpreted 
as a relativistic generalization of purely rotational fiuid fiow. Here, they appear 
exclusively associated to pure shear motion. The vorticity for this mode is zero, 
since its relevant covariant velocity component U2 = (indeed, from Euler equations: 
U2 = -u^ + ho2 = 0). 

We can now discuss the two main aspects of this problem: the first aspect con- 
cerns the physics of the instability, and the second the generation and propagation 
of stable waves. 



Starting with the first aspect, one notices from equation ([1^ , that the instability 
sets in through an a; = mode, with a critical wave number kc = a/SttGpo/c^, where 
the speed of light c has been restored. In terms of wavelength, above the critical 
value Ac = ^Tr/kc, i.e., for 

/ 2 \ 1/2 

the system is unstable, A = Ac being the marginal case. This instability tends to 
shear in the material along the node of the perturbation. Indeed, under the effect of 
the wave the particles of the fluid get an acceleration r2oo = ^02 that produces the 
velocity = /io2- As shown in Figure ^ one sees that part of the material moves to 
one side, part to the other, whereas it stays at rest at the node where the velocity 
in the y— direction, u"^, is zero. 

A very important point which should be made explicitly clear is that at such 
critical wavelengths there are equilibrium configurations which are not at all dy- 
namically unstable. Indeed, in a previous paragraph we have put the following 
limits on the equilibrium configuration, i -C Csl\fGp and iL <^ c^/Gp. To have 
L > Ac ~ cj \JGp while the above conditions are still satisfied, I has to satisfy both 
limits £ <C CgXc/c and ^ <C A^/L. Thus, for sufficiently small thickness there are 
well defined equilibria. 

In order to better understand such an instability it is advisable to rearrange 
equation ( p!7|) above and write it as /z > vrm^/A, where fj, = aX = poiX is the mass 
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Figure 1: A helicity-1 shear wave of gravity travels in the 2— direction through a 
plane parallel slab and induces a motion on the fluid in the ?/— direction. The static 
gravitational potential is a function of the x coordinate. 

per unit length in the direction and = c^i/2G is the reference relativistic mass 
associated with the thickness length i. Thus m^/A can be seen as the relativistic 
reference mass per unit length in the z— direction. The instability criterion (|T^ 
says that when the mass per unit length in the direction of the perturbation is 
about greater than the reference mass per unit length in the z— direction then the 
momentum in the |/— direction is enough to overcome the pull back from the mass 
along z and the material breaks apart. 

This instability is of course different from Jeans instability but not quite in 
what regards the gravitational field dynamics. Equation ([12|) can be written as 
n^02 = StcGST^"^, because = hQ2, and if one performs a Lorentz boost along the 

direction, it is equivalent to Oh'^Q = —SnG 6T'^^ in the new frame, which is com- 
parable to the equation governing the response of the gravitational field to density 
fluctuations. What is different is that they are not accompanied by perturbations 
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in the fluid proper density and pressure as in the Jeans, acoustic case, but rather 
by transverse displacements of the fluid. For such a plane layer and density wave 
perturbations of the form exp{—iujt + ikz), Jeans instability sets in when k < kj 
with kj = ^yAnGpo/cl, or X > Xj with Xj = a/vtc^V^Po- The material collapses 
into a set of plane-parallel slabs distributed perpendicularly to the 2;— direction and 



uniformly along the z— direction p3|, 



The instability we have found here sets in through the interaction of shear waves 
and matter. A helicity-1 shear wave travelling through the plane parallel layer in 
the 2— direction with A > Ac will induce an instability in the matter and in the 
gravitational fleld itself. For a planar slab characterized by a surface density a and 
height i one can write the critical wavelength as Ac = \/ TTc^i/2G(7. One sees that 
Ac can be relatively small for high enough surface densities or very thin systems. 
Conversely, if such instabilities are not detectable within the system, one can put a 
lower limit on the ratio i/a; for a slab with size L these instabilities do not show 
up if Ac > L, i.e., £/a > (2G/7rc^)L^. An estimate of Ac for our Galaxy yields 
Ac — 16 000Kpc showing that it is not thin enough to present such instabilities. 
However, for extended and very thin systems this instability is applicable and could 
set in. 

Within plane parallel symmetric conflgurations, other unstable modes should 
occur, although we have been able to solve the system consistently only for the 
ho2 — mode because we have restricted this study to plane waves. It would be 
interesting to find a consistent solution for the mode involving /iqo- 

We have considered here plane parallel symmetry but other type of configura- 
tions, such as disk-like, should also possess such instabilities. For disk-like con- 
figurations some other features and new global instabilities could set in. Jeans 
instabilities show surprises in the disk-like case, since for high enough values of the 



sound speed a disk is stable to all axisymmetric perturbations [25, Eq], a feature 



that does not appear in the spherical and planar cases. Domain walls and cosmic 
strings should be considered. In a preliminary calculation, we have found that do- 
main walls also present this type of instability to wave perturbations of the gravity 
and scalar field [^. On the other hand, the spherical case is less interesting since 
here Ac ~ R^^'^/rI^'^, where Rs is the Schwarzschild radius of an object of radius 
R. Thus, the instability sets in when i? ~ i?^ ~ Ac, and other stronger dynamical 
instabilities must already be present. 

We now discuss the second aspect, i.e., the generation and propagation of stable 
waves of helicity-1. In order to estimate the energy fiux of these waves we use the 
expression for the effective stress-energy t^i, tensor given in |18] and apply to our 
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case. It yields 



iau — , „ ^ < —^02,(4^02,1/ + ^23.u^23,i/ >— ""TTPT ^ ^02,(t^02,i/ > , (18) 
iovrG ' Zfc^ 

where the symbol < > denotes an average over one wave period. Defining = 
\u\ exp{—iut + ikz), with \u\ being the amplitude of the velocity in the direction, 
we find 

= IpoH'v , (19) 

where v is the velocity of propagation of the wave. It is now difficult to estimate at 
this point a realistic energy fiux, given one has to deal with the propagation as well 
as the generation of such a wave. One possibility is to have a source of mass M in 
an oscillatory (orbital) movement of radius R, inducing tidal motions on particles 
at a distance r. This induces a tidal velocity to the particles which we identify with 
\u\, given then by \u\/c ~ y/RgR/r, where Rg = 2GM/c^. If we put Rg ~ 10~^R, 
R ~ 10~^r and po = Ig/cm'^ we obtain \u\/c ~ 10~^ and, for f ~ c, an energy 
fiux t^^ ~ 10^° erg/(s cm^). Perhaps, the most important point is to notice that 
this tidal perturbation propagates away through the matter at typically the speed 
of light according to the dispersion relation (|T^ and properties we have derived for 
this helicity-1 wave. One can then think that randomly distributed waves of the type 
considered, should induce thermal motions in the matter, which would be detected 
as velocity dispersions in the stars or particles of the system. These velocities, in 
turn, would give the impression of the existence of an extra Newtonian gravitational 
potential which, for instance, could be interpreted as coming from dark matter, or 
from modifications of the gravitational potential itself. 

The propagation of these helicity-1 shear waves and the instabilities themselves 
we have been studying here can be considered wave detectors in potential and deserve 
close attention. 
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